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(Dated: February 7, 2008) 

Using events in which one of two neutraJ-B mesons from the decay of am T{4S) resonance is fully 
reconstructed, we set limits on the difference between the decay rates of the two neutral-B mass 

cigcnstates and on CP, T, and CPT violation in B^B'^ mixing. The reconstructed decays, compris- 
ing both CP and flavor eigenstates, are obtained from 88 million T(45) BB decays collected 
with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. Wc determine six 
independent parameters governing oscillations (Am, AF/F), CPT and CP violation (Rez, Imz), 
and CP and T violation (ImAcp, where Xcp characterizes B° and B° decays to states of 

charmonium plus Kg or ■ The results are 

sgn(ReAcy)AF/r = -0.008±0.037(stat.) ± 0.018(syst.) [-0.084,0.068], 

\q/p\ = 1.029±0.013(stat.) ±0.011(syst.) [ 1.001,1.057], 

(ReAcy/|Acp|)Rez = 0.014±0.035(stat.) ± 0.034(syst.) [-0.072,0.101], 

Imz = 0.038±0.029(stat.) ±0.025(syst.) [-0.028,0.104]. 

The values inside square brackets indicate the 90% confidence-level intervals. The values of Im \cp 
and Am are consistent with previous analyses and are used as cross-checks. These measurements 
are in agreement with Standard Model expectations. 

PACS numbers: 13.25.Hw, 12.15.Hh, 14.40.Nd, ll.30.Er 



I. INTRODUCTION AND ANALYSIS 
OVERVIEW 

The mass difference Am between the mass eigen- 
states has been measured with high precision at B- 
factory experiments [1 4], and CP violation has been 
observed in neutral- i?-meson decays to states hke 
J/tpK^ [5, 6]. However, our knowledge of other aspects 
of neutral- B-meson oscillations is meager. In this paper, 
wc provide direct limits on the total decay-rate difference 
AF between the B° mass eigenstates, and on CP, T, and 
CPT violation due to oscillations alone. 

In the Standard Model, the ratio AT /Am is of or- 
der ml /ml and thus quite small. Recent calculations 
of Ar/r, including 1/mb contributions and part of the 
next-to-leading order QCD corrections [7, 8], find values 
in the approximate range —0.2% to —0.3%. Existing lim- 
its for I AF/F] [9, 10] are relatively weak (~ 20%). The 
large data sets available at asymmetric-energy B facto- 
ries provide an opportunity to look for deviations from 
the Standard Model. 

The CP-violating asymmetry observed in neutral-S- 
meson decays to states like J/ipK^ is due to the in- 
terference between decay amplitudes to a CP eigen- 
state with and without mixing. CP violation in mixing 
alone leads to different rates for the transitions B^ B^ 
and -B" B'^ . This can be measiued, for example, by 
comparing the decay rates to £^£^X and from 
semileptonic decays of pairs of neutral-B mesons arising 
from the T(4S') [11]. The only semileptonic decays gen- 
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erated by first-order weak interactions are B^ i'^vX 
and B'^ — > £~VX and the CP invariance of strong and 
electromagnetic interactions guarantees that these have 
equal rates. As a result, any asymmetry in the dilepton 
rates can be ascribed to CP violation in mixing. While 
CP violation in mixing is suppressed in the Standard 
Model [8, 12, 13], additional virtual contributions from 
new physics could obviate this suppression. Similarly, 
new physics may introduce additional intrinsic T viola- 
tion or even CPT violation in mixing. It is these possi- 
bilities for the breaking of discrete symmetries in mixing 
itself that we address in this analysis using nonleptonic 
decays that are completely reconstructed. 

The behavior of neutral-B mesons is sensitive to CPT 
violation [14-16]. A theorem [17] founded on general 
principles of relativistic quantum field theory states that 
the CPT symmetry holds for any local field theory satis- 
fying Lorentz invariance. The CPT symmetry is the only 
combination of C, P, and T that is not known to be vi- 
olated. Nevertheless, it is possible that CPT symmetry 
could fail at short distances [18]. Strict constraints on 
CPT violation have been obtained in the neutral-kaon 
system [19]. Limits in the B-meson system have been 
obtained previously [4, 20]. 

To measure AF and CP, T, or CPT violation, we ob- 
serve the time dependence of decays of neutral- i? mesons 
produced in pairs at the T(45) resonance. The usual ap- 
proach to mixing and CP analyses [1-6] allows for expo- 
nential decay, modulated by oscillatory terms with fre- 
quency Am. These analyses neglect the difference AF 
betwc!c;n the decay rates of the two mass eigenstates, 
which would introduce terms with a new time dependence 
exp(±AFt/2). Violation of CP, T, or CPT in the mixing 
of the neutral-_B mesons would modify the coefficients 
of the various terms involving exponential and oscilla- 
tory behavior. To detect these potential subtle changes 
requires precision measurements of the decays, detailed 
consideration of systematic effects, and thorough treat- 
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ment of coherent production of neutral- B-meson pairs 
from the 7(45). 

This analysis is based on a total of about 88 million 
r(4S') BB decays collected with the BABAR detector 
at the PEP-II asymmetric-energy B Factory at the Stan- 
ford Linear Accelerator Center. There, 9.0- GeV elec- 
trons and 3.1- GeV positrons annihilate to produce the 
BB pairs moving along the e~ beam direction (z-axis) 
with a Lorentz boost of l3"f w 0.55. This boost makes 
it possible to measure the proper-time difference At be- 
tween the two B decays. We fully reconstruct one me- 
son from its decay to a flavor eigenstate (-Bfiav) or to a 
CP eigenstate {Bcp) composed of charmonium and ei- 
ther a or K^. We denote the flavor and CP cigen- 
states jointly by Brec- The remaining charged particles 
in the event, which originate from the other B meson 
(^tag), are used to identify ("tag") its flavor as B° or 
B^. Not all events can be tagged, but the untagged 
events are also used in the analysis. The time diff'er- 
ence At = t-^cc — ttag ~ A.z/{/3"fc) is determined from 
the separation Az along the boost direction of the decay 
vertices for the fully reconstructed B candidate and the 
tagging B. 

A maximum-likelihood fit to the time distributions of 
tagged and untagged, flavor and CP eigenstates deter- 
mines six independent parameters (sec Sec. II) govern- 
ing oscillations (Am, AT/T), CPT and CP violation 
(Rez, Imz), and CP and T violation (ImAcp, \q/p\), 
where \cp is the usual variable used to characterize the 
decays of neutral-S mesons into flnal states of charmo- 
nium and a K'^ or K^. The values of ImA^p and Am 
are used as cross-checks with the earlier BABAR sin 2/3 
result [5], obtained with the same dataset, and with pre- 
vious B- factory measurements of Am [1-4]. All the pa- 
rameters are explicitly defined in Sec. II. 

The analysis presents several challenges. First, the res- 
olution for At is comparable to the B lifetime and is 
asymmetric in At. This asymmetry must be well under- 
stood lest it be mistaken for a fundamental asymmetry 
we seek to measure. Second, tagging assigns flavor in- 
correctly some fraction of the time. Third, interference 
between weak decays favored by the Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix and those doubly- 

I 



WfJ-,L = 'mH,L - -^^H,L = m - 



where the real part of the square root is taken to be 
positive and where we define 

m=i(Mn+M22) , T = i(rii + Fas) , 

5m = Mil - M22 , 6r = Fii - F22 . (2) 



Cabibbo-suppresscd (DCS) cannot be neglected. Fourth, 
direct CP violation in the Bcp sample could mimic CP 
violation in mixing and must be parameterized appropri- 
ately. Finally, we have to account for possible asymme- 
tries induced by the differing response of the detector to 
positively and negatively charged particles. In resolving 
all of the above issues we rely mainly on data. 

This paper provides a detailed description of the anal- 
ysis published in Ref. [21], and is organized as follows. 
In Sec. II we present a general formulation of the time- 
dependent decay rates of B°B° pairs produced at the 
T{4:S) resonance, including effects from the decay-rate 
difference, possible CP and CPT violation in mixing, 
and interference effects induced by DCS decays. We de- 
rive the expressions for B decays to flavor and CP eigen- 
states. In Sec. Ill we describe the BABAR detector. After 
discussing the data sample in Sec. IV, we describe the B- 
flavor tagging algorithm in Sec. V. Sec. VI is devoted to 
the description of the measurement of Az and to the de- 
termination of At and its resolution function. In Sec. VII 
we describe our log-likelihood function and the assump- 
tions made in the fit. The results of the fit are given 
in Sec. VIII. Cross-checks arc discussed in Sec. IX and 
systematic uncertainties are presented in Sec. X. The 
results of the analysis are summarized and discussed in 
Sec. XI. 



II. GENERAL TIME-DEPENDENT DECAY 
RATES FROM r(4S) B°B° 

The neutral- B-meson system can be described by the 

effective Hamiltonian H — M — ir/2, where M and T are 
two-by-two Hermitian matrices describing, respectively, 
the mass and decay-rate components. CP or CPT sym- 
metry imposes that Mn = M22 and Fn = F22, the in- 
dex 1 indicating B" and 2 indicating B°. In the limit of 
CP or T invariance, F12/M12 = F21/M21 = T^^/M^^, so 
F12/M12 is real. These conditions do not depend on the 
phase conventions chosen for the B° and B°. The masses 
'nH,L and decay rates Th,l of the two eigenstates of H 
form the complex eigenvalues coh,l 




I 

Assuming CPT invariance {Sm = 0, (JF = 0), and antici- 
pating that |AF| <^ Am, we have 



Am = mH - rriL ~ 2|Mi2| , 
AF = Th-Tl^ 2IM12I Re(Fi2/Mi2) . (3) 
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Here wc have taken Am to be the mass of the heavier 
eigenstate minus the mass of the hghter one. Thus, AF 
is the decay rate of the heavier state minus the decay rate 
of the hghter one and its sign is not known a priori. 

With CPT symmetry, the hght and heavy mass eigen- 
states of the neutral- B-meson system can be written 



\Bl) = p\B^)+q\B^) 
\Bh) = p\B°) - 



where 



g 
P 



M- 



12 



The magnitude of q/p is very nearly unity: 

2 n 



Im- 



12 



M- 



12 



(4) 



(5) 



(6) 



In the Standard Model, the CP- and T-violating 

quantity jg/pp — 1 is small not just because |ri2| is 
small, but additionally because the CP-violating quan- 
tity Im(ri2/Mi2) is suppressed by an additional factor 
{ml — ni'^)/ml w 0.1 relative to |ri2/Mi2|. Violation 
of CP is not possible if two of the quark masses (for 
quarks of the same charge) are identical, for then we 
could redefine two new quark states with equal masses 
so that one of them did not mix with the two remaining 
states. The mixing among two generations would be in- 
adequate to support CP violation. When the remaining 
Standard Model factors are included, the expectation is 
|Im(ri2/Mi2)| < 10-3 [8, 12, 13]. 

CPT violation in mixing can be described conveniently 
by the phase-convention-independent quantity 

Sm-^ST 



Ml 



■ ri2) (M- 



12 



2 ^ 12) 



m - 



dm 



ST 



Am - ^Ar 



(7) 



The generalization of the eigenstates in Eq. (4) when we 
account for CPT violation can be written 



\Bl) = pVT^lB"" 
\Bh) 



qVTT^lB^) 
pVT+^\B°) - qVT^\B°) 



(8) 



where we maintain the definition of q/p given in Eq. (5). 
The result, when time evolution is included, is that states 
that begin as purely B° or B^ after a time t will be 
mixtures 

|B°phys(t)) = [9+{t)+zg_{t)] \B°) 

-v/r^%_(t)|BO) , (9) 

I^VysW) = [9+{t)-Z9.{t)]\B°) 

-Vl^^g.{t)\B^) , (10) 



where we have introduced 



Invariance under CP or under T requires that 

|(i?0|B"phys(0)l = l(S°|B%hys(i))| ; 



(11) 



(12) 



i.e., \q/p\ = 1, which is guaranteed by Im(ri2/Mi2) = 0. 
Table I shows the constraints on \q/p\ and z for the dif- 
ferent possible symmetry scenarios. The Standard Model 
corresponds to the second configuration (CPT symme- 
try, with CP and T violated). Note that two of these 
scenarios are degenerate. With CP symmetry in B'^B'^ 
oscillations, this experiment cannot distinguish between 
T and CPT both being conserved or violated. 



TABLE I: Constraints on \q/p\ and z due to CP, T, and CPT 
symmetries in B°W oscillations. 





CPT 


CPT 




CP, T 




gp , T 


CP, jr 






= 1 




= 1 


= 1 


^1 


z 


= 


= 


^0 


= 


^0 



A. Effects of Coherence 

At the T{4:S) resonance, neutral-B mesons are pro- 
duced in coherent p-wave pairs. If we subsequently ob- 
serve one B-meson decay to the state /i at time to = 
and the other decay to the state /2 at some later time 
t, we cannot in general know whether /i came from the 
decay of a B° or a and similarly for the state /2. If 
Ai^2 and Ai^2 are the amplitudes for the decay of i?" and 
B°, respectively, to the states /i and /2, then the overall 
amplitude is given by 



= o,+g+{t) + a-g-(t) , 

where 

a+ = -A1A2+A1A2 , 
a_ = Vl-z2 



(13) 



^AiA2 - -A1A2 



+ z [A1A2 + A1A2] . 

(14) 



Using the relations 
1 



9±{t)? = ^ e-" [cosh(Ari/2) ± cos(Amt)] (15) 



and 



gl{t) 9-{t) = --e-'* [sinh(Art/2) + i sin(Amt)] , 



(16) 
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we find the decay rate 

-T^ oc e"^l*li ic+cosh(Art/2) + ic-cos(Amt) 

Q.i' I Zi Zi 

-Res sinh(Ari/2) +Ims sin(ATOi) \ , (17) 



where 



c± 



(18) 



The absohitc vahic in the leading exponential in Eq. (17) 
is introduced for later convenience. 

Now let us take /i = /tag to be the state that is in- 
completely reconstructed and that provides the tagging 
decay, and /2 = /rcc to be the fully reconstructed state 
(flavor or CP eigenstate). Then we have t = ti. 
and Eq. (14) becomes 



t 



tag 



tag-'^rec ~t~ -^tag-'^rec j 



a_ = Vl-z^ 

q " p 

+Z I^Ajag^rec + ^tag^rec] • 



Pa, a -^a. a 

-^t ag -^rec ^ ^tag -^rec 



(19) 



If instead the tagged decay occurs second, we would 
need to redefine t, and a_ by interchanging the labels 
"tag" and "rec" . This would amount to the replacements 
t —t, a+ —a+, and a_ ^ o_. However, we see that 
Eq. (17) is actually unaffected by these changes and that 
we can instead retain the definitions t — t,cc — hag and 
those of Eq. (19). Thus, Eqs. (17)-(19) apply indepen- 
dent of the order of the decays of the tagged and fully 
reconstructed B mesons. 

A fully reconstructed flavor state cannot always be un- 
ambiguously associated with either B° or B°. DCS de- 
cays, such as B^ — > £)+7r~, occur at a rate suppressed 
by roughly K.VJV^.V^J^ « (0.02)^. Although this 
can be neglected, interference between favored and sup- 
pressed amplitudes is reduced by only a factor of approx- 
imately 0.02 [22], and must be taken into account. 

Tagging cannot l)e done perfectly, largely because the 
tagging state is incompletely reconstructed. We account 
for this by measuring the wrong-tag probability from the 
data. However, even if our tagging were perfect in prin- 
ciple, it would be afflicted with the same complication 
from DCS decays as the fully reconstructed state. The 
full expressions for the real coefficients c± and the com- 
plex coefficient s, containing the DCS amplitudes, are 



c± = 



^rec^tag ^rec^tae I i l^l l^rec^tae + ^rec^tag I i|l ^ I 



Lrec-^tag -^rec-'^tag 

±2Re 



Ia a -^A X 

-^rec-^tag -^rec-^tag 



\/l — ^"^rcc^tag "^rcc^tag^ (^rec^tag + ^rec^tag) I" 



i^-'^-tag -^rec-^tag) 



^^-^rec-^tag ^-^rec-^tag^ H" ^ (-^rec-^tag H~ -^rec-^tag) ^ • 



(20) 
(21) 



TABLE II: The coefficient c+ from Eq. (20), evaluated to leading order in the small quantities z. As/, ^Bt, ^gf, and A^j. If 
the tagging state for a B" tag is /st, then the tagging state for a B" is the CP-conjugate state, fg^, and siniilaiiy for the fully 
reconstructed states. The decay amplitudes are Aat = {fBt\H\B°), Aat = (/st|7i|B°), Ag^ = (/bJ'W|-B°), Ag^ = {} bi,\'H\B°) , 
and similarly for rec = Bf, Bf, CP. 



Btag 




c+ 


B° 


5° 


\ABt 




B° 


B" 


\ABt 




B° 


B" 


\ABt 






B" 




'\ABf?Vl/P? 




Bcp 


\ABt 


''\Acp\' \p/q\'[l + |Agp|^ -4ReA(y ReAist -F2RezReAcy -2ImzImAa'] 


B° 


BcP 


Pa 


'^\Acp\^ [1 + \\cp\^ -4ReAcy ReAgj - 2RezReAcp - 2ImzImAcp] 



Terms proportional to A^ec^tag and Arec^tag are as- sociated with decays with no net oscillation between 
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TABLE III: The coefficient c- from Eq. (20), evaluated to leading order in the small quantities z, As/, Ast, Agy-, and Agj. See 
caption of Table II for the definition of the various quantities. 



-'-'tag 




c_ 


B° 


B° 


-\ABt\MBf?\p/q? 


B° 


B° 


\ABt?\Asf\^ 


B° 


B° 


\Am?M 


B° 


B° 


-\AEt?\AEj?\llP? 


B° 


BcP 


\ABt\^\Acp\^ \p/q\^[-l + [AgpI^ -4ImAc2pImAst - 2RezRe Acy + 2ImzImAcp] 


B° 


BcP 


Pstl^l^GPl^ [1 - |Acp|^ +4ImAcp ImAs^ + 2RezReAcy + 2ImzImAc2p] 



TABLE IV: The complex coefficient s from Eq. (21), evaluated to leading order in the small quantities z, Ab/, Ab*, ^gj, and 
Xgf. See caption of Table II for the definition of the various quantities. 



Stag 




s 




B" 


\ABt 


Mbi 


^\p/q\^[\*B^-\*Bf] 






B" 


\ABt 


'\Abs 


''[\Bt-\Bi-A 






B" 


lAst 




'pBt-AB/+z] 






B" 


lAgt 


''\ABf 






B« 


Bcp 


\ABt 




''\plq\'[\\cp\'\Bt-\lT 


+ A|5i - IAcpI^'z] 


B° 


BcP 


lAst 









the two neutral-i? decays, while terms proportional to 
(gr/p)ArecAag and (p/5)j4rec^tag represent a net oscilla- 
tion. 

We characterize each final state / through the param- 
eter 

where / can be "rcc" (which can be itself "flav" or "CP" ) 
or "tag". In the absence of DCS decays, Aflav (i-e., Arec 
when the reconstructed state is a flavor eigenstate, not a 
CP eigenstate) and Atag would be either zero or infinite. 
With a contribution from DCS decays they arc non-zero 
and finite. 

If the reconstructed flavor state /flav is ostensibly a 
(hereafter indicated as Bf to avoid ambiguities with the 
tag state) then lAiod = \^Bf\ ^ 1- Conversely, if the re- 
constructed state appears to come from a B^ (indicated 
as Bf), then |Ajjj| ^ 1, and it is convenient to intro- 
duce = 1/A^^. The pattern for the tagging state 
("tag"= Bt, Bt) is similar. If the reconstructed state is a 
CP eigenstate, then |Arec| = |Acp| is of order unity. 

In practice, terms quadratic in z or in a small A/ are 
not important. The expressions for c± and s when only 
linear terms in small quantities are retained are shown 



in Tables II, III, and IV. The analysis uses the full 
expressions, without simpliflcation. 

It is appropriate to assume that the decays to flavor 
eigenstates we consider arc dominated by a single weak 
mechanism: h cud. While we can flnd a mechanism for 
b cud (which is a DCS process), there are no alterna- 
tive first-order weak processes that produce cud from a b 
quark. Then even if there are several contributions to the 
decay, each possibly with its own strong phase, the CP- 
conjugate decay differs only by changing a single common 
weak phase so that = \Aqj:\, = l^g/l (and 

similarly for tagging states). In fact, even if this assump- 
tion is not rigorously true, any violation will be absorbed 
in tagging and reconstruction efficiencies, which are de- 
termined from the data, as described in Sec. VII. These 
equalities relate the four permutations that arise from 
the tag and reconstructed state being either B° or B°. 



B. Ensembles of States 

In principle, every hadronic final state fh has a differ- 
ent Aft, which can be written as A^ = |A;j|e~"^^'^''°e~"^^ , 
where <^^™" and (p'jf'^ are strong (CP-even) and weak 
(CP-odd) phases that arise from the ratio of the ampli- 
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tudes of the and B'^ decays to fh- Assuming that 
there is a single weak phase involved, the CP-conjugate 
state will have \ = b/gplAf^le-^^^^^'e'^S''. 

If we sum squares of amplitudes over a collection F of 
flavor states that are ostensibly B^, the terms that do 
and do not contain Aflav are of the form 

^ lA^l'Xh and Yl 1^"!' ' (23) 
fheF fheF 



so we can define an effective As/ by 
Similarly, for flavor states that are ostensibly 



(24) 



hi '^h 



(25) 



The two complex numbers A_b/ and A^j encapsulate the 
effects due to DCS decays in the fully reconstructed B 
decay, as long as the terms quadratic in A/j and A^, sup- 
pressed by roughly |V;*(,Kd/K6K.dP ~ (0-02)', are omit- 
ted. 

The same argument applies to tagging states. If the 
collection of states contributing to a B^ or B^ tag is a 
then 



Asf — 



^Bt — 



A. 



(26) 
(27) 



In practice, we do not use separate Atag parameters for 
each tagging category a (i.e., each collection of states 
of similar character, as described in Sec. V), but simply 
one for B° and one for B^, setting aside the lepton tag 
category, which is free of DCS decays. This treatment 
is flexible enough to incorporate the DCS-decay effects 
that can mimic the asymmetries we seek in the analysis. 

Henceforth, expressions like Aflav and Atag refer to an 
appropriate sum over observed states. The summation 
over states fh in a tagging category should be thought 
of as extending over those states that are reconstructed 
as belonging to the given category. In this way, we in- 
corporate implicitly the tagging efficiency of each state 
fh- The reconstruction efficiency is incorporated in an 
analogous fashion into Aflav 

Data from directly related CP final states like J/ipK^, 
with r]cp = — 1, and J/tpK^, with rjcp = +1, where r]cp 
is the CP eigenvalue of the final state, can be combined 
by assuming that their time distributions arc identical, 
except for the factor rjcp. We use a single parameter 
Acp obtained multiplying Eq. (22) by r/cp- We assume 
rcp = \Acp/Acp\ = 1 as expected theoretically at the 
lO"*^ level [23] and as supported experimentally by: 



i) the average of _B-factory measurements of states 
of charmonium and or K^, from which it has 
been obtained rcp = 0.949 ±0.045 [5, 6], when AT, 
\q/p\ — 1 and z are assumed to be zero; 

ii) the average of CLEO and BABAR measurements 
of the CP asymmetry in the charged mode 
B^ J/ipK^ , from which it is found r(jp_j/^j(± = 
1.008 ± 0.025 [24, 25], combined with isospin sym- 
metry to relate with the CP final states [26]. 



C. Sensitivity of Distributions to Parameters 

From Eq. (17) and Tables II, III, and IV, it can be 
seen that while ImAcp, Imz, \q/p\, and rcp are unam- 
biguously determined, Rez appears only in the prod- 
uct Re Acp Re z or else is suppressed by the small fac- 
tor Ar/r. Similarly, the sign of AF cannot be deter- 
mined separately from the sign of ReAcp since AF al- 
ways appears multiplied by ReAcp in its dominant con- 
tribution. Its value is known only through Re Xcp = 
±\/|Acp|^ — (ImAcp)^, where the choice of sign could 
be made by a separate measurement that directly deter- 
mines the sign of Re Xcp- As a result, the parameters that 
can be determined by this analysis are sgn(Re Acp)AF/F, 
\q/p\, (Re Acp /I Acp I) Rez, Imz, ImAcp/jAcpj, rcp, Am, 
and F. In practice, we fix rcF and F in the nominal fit, 
and vary them for systematic studies. 

Data for final states that are CP eigenstates and 
those that are fiavor eigenstates are both needed for 
the analysis, as shown in Table V. The sensitivity to 
(Re Acp /I Acp I) Rez and ImAcp/jAcpj is provided by the 
decays to CP eigenstates Bcp, for which the accompa- 
nying t dependence is even for the former and odd for 
the latter. The -Bflav sample contributes marginally to 
these parameters because it lacks explicit dependence on 
Im Acp /I Acp I and the dependence on Rez is scaled by 
the sinh(AFt/2) term, which is small for small AF. 

In contrast, the parameters \q/p\ and Imz (and Am) 
are determined by the large -Bflav sample, where the for- 
mer is associated with a t-even distribution and the latter 
with a t-odd distribution. For small values of AF/F, the 
determination of AF/F is dominated by the Bcp sam- 
ple, despite the smallness of this sample compared to the 
-Bflav sample. This is because in the flavor sample the 
leading dependence on AF is proportional to AF^ , while 
in the CP sample it is proportional to AF. The contribu- 
tion of sinh(AFt/2) is the same for both 5° and 5° tags, 
so events that cannot be tagged may be included in the 
analysis to improve sensitivity. The Bcp sample is also 
sensitive to the sign of AF/F (up to the sign ambiguity 
from Re Acp). 

Overall, the combined use of the Bflav and Bcp samples 
provides sensitivity to the full set of physical parameters, 
since they are determined either from different samples, 
or from different t dependences. 

As we show in Tables II, III, and IV, if the recon- 
structed state is a flavor eigenstate, the DCS-decay effects 
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TABLE V: Dominant dependence of the time distributions 
on the physical parameters measured with fully reconstructed 
flavor and CP states. Sensitivity is specific to terms in the 

time dependence that are either t-evon or t-odd. The flavor 
sample is much larger than the CP sample. 



Parameter 


B 

t-even 


t-odd 


B 

t-even 


CP 

t-odd 


\qIp\ 

Am 


X 
X 








Imz 




X 






(ReAGp/|Acp|)Rez 
rcp 






X 
X 




sgn(ReAGp)Ar/r 
Im Acy/ Acy 








X 
X 



in tagging are negligible except in the sin(Amf) term, for 
the other terms are suppressed by both a power of Aflav 
and a power of Atag- Conversely, if the reconstructed 

state is a CP cigcnstatc with \Xcp\ ~ Ij the effects from 
DCS decays are confined to the terms even in t. 



III. THE BABAR DETECTOR 

The BABAR detector is described in detail else- 
where [27], so here we give only a brief description of 
the apparatus. 

Surrounding the beam-pipe is a five-layer silicon ver- 
tex tracker (SVT), which gives precisely measured points 
along the trajectories of charged particles as they leave 
the interaction region. Outside the SVT is a 40-layer drift 
chamber (DCH) filled with an 80:20 helium-isobutane 
gas mixture, chosen to minimize multiple scattering. 
Charged-particle tracking and the determination of mo- 
menta through track curvature rely on the DCH and SVT 
measurements in the 1.5-T magnetic field generated by a 
superconducting solenoid. The DCH and SVT measure- 
ments of dE/dx energy loss also contribute to charged- 
particle identification. 

Surrounding the drift chamber is a novel detector of 
internally reflected Cerenkov radiation (DIRC), giving 
charged-particle identification in the central region of 
the detector. Outside the DIRC is a highly segmented 
electromagnetic calorimeter (EMC) composed of CsI(Tl) 
crystals. The EMC is used to detect photons and neu- 
tral hadrons through shower shapes and is also used to 
identify electrons. Finally, the flux return of the super- 
conducting coil surrounding the EMC is instrumented 
with resistive plate chambers interspersed with iron for 
the identification of muons and neutral hadrons (IFR). 

A detailed Monte Carlo program based on the 
GEANT4 [28] software package is used to simulate the 
BABAR detector response and performance. 



IV. DATA SAMPLES AND B-MESON 
RECONSTRUCTION 

From a sample of about 88 million T(4S') BB de- 
cays, we select events in which one of the B mesons is 
completely reconstructed in either a neutral or a charged 
hadronic final state, using the same criteria used for the 
BABAR sin 2/3 measurement [5] and for measiirements of 
Am using hadronic final states [1]. Neutral- i? mesons 
are reconstructed in either a flavor (Bflav) or CP (Bcp) 
eigenstate. The charged- B- meson decays are used as con- 
trol samples in the cross-checks described in Sec. IX B. 
The decay modes used for the flavor sample, the CP sam- 
ple, the control samples are displayed in Table VI. De- 
tails on charged particle and neutral reconstruction, par- 
ticle identification and reconstruction of B mesons can 
be found in Sees. II and HI in Ref. [29]. 

TABLE VL The flavor, CP, and control sample decay modes 
used in this analysis. The J/ip is always identified in the e"'"e~ 
or 11^ modes. The is reconstructed only in tt+tt+tt". 
The K% is identified in the 7r"'"7r~ mode, except otherwise 
specified. All charge-conjugate decay modes are included im- 
plicitly. 



Samples 


Decay modes 


-Sflav 


B° ^ D—i<:+{p+,a+) 

D*- -^T}\- 

B° D-7v+{p+,a+) 

D- K+T:-it-,K°s'K- 
K''° K+TT- 


Bcp 


B°^J/^K°s 
B° ip{2S)K°s 

ip{2S) e+e",/i+/i", 
Xci J/tp 7 

B°^J/i;K° 


Control 


B+^D^*'\+ 

B+ J/tljK+ 
B+ -> ^p(2S)K+ 
B+ ^ XciK+ 
B+ -> J/ipK*+ 

K*+ K°sn+ 



We select Bflav and Bcp candidates by requiring that 
the difference AE between their energy and the beam 
energy in the center-of-mass frame be less than Sa from 
zero, where a is the resolution on AE. The AE' resolu- 
tion ranges between 10 and 50 MeV depending on the 
decay mode. For Bflav modes and Bcp modes involving 
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{Bqpj^o ) , the beam-energy substituted mass must be 

greater than 5.2GeV/c^. The beam-energy substituted 
mass is given by 



mEs 



'(s/2 + p, •pb)2 



Pi 



(28) 



where s is the square of the center-of-mass energy, Ei and 
Pi are the total energy and the three-momentum of the 
initial state in the laboratory frame, and ps is the three- 
momentum of the B candidate in the same frame. In 
the case of decays to J/ijj {Bqpj^o ) , the direction 
is measured but its momentum is only inferred by con- 
straining the mass of the J/ip candidate to the known 
B^ mass. As a consequence, there is only one parame- 
ter left to define the signal region, which is taken to be 
\AE\ < lOMeV. 

Fig. 1 shows the toes distribution for the -Bflav and 
Bqpj^o samples and the AE distribution for the Bfjp^o 
candidates, before the vertex requirements (see Sec. VI). 
The combinatorial background in the toes distribu- 
tions is described by the empirical ARGUS phase-space 
model [30] and the signal by a Gaussian distribution. 
The combinatorial background consists of random com- 
binations of tracks from continuum and BB sources. 
The former events are dominantly "prompt" , that is, 
the observed particles point back to the interaction 
point, whereas the latter events are dominantly "non- 
prompt" , with particles pointing back to separated ver- 
tices. Charmed particles, either from continuum or from 
B-meson decays, contribute to non-prompt background. 
A small background due to other B decays (not shown 
in Fig. 1) also peaks at the B mass. The background in 
the J/ipKl channel receives contributions from other B 
decays with real J/ijj mesons in the final state, and from 
events with fake J/^ mesons constructed from unassoci- 
ated leptons or from misidentified particles. 

After completely reconstructing one B meson, the rest 
of the event is analyzed to identify the flavor of the op- 
posite B meson and to reconstruct its decay point, as 
described in Sees. V and VI. 

Using exactly the same requirements, we analyze 
GEANT4-simulated samples to check for any biases in 
the event selection and extracted parameters. The Monte 
Carlo samples are also used in studies of detector re- 
sponse and to estimate some background sources. The 
values of the B oscillation and CP-, T-, and CPT- vio- 
lating parameters assumed in the simulations are similar 
to those measured in the data. We use additional samples 
with significantly different values to check the reliability 
of the analysis in other regions of the parameter space. 



V. FLAVOR TAGGING 

The tracks that are not part of the fully reconstructed 
B meson are used to determine whether the -Btag was 
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FIG. 1: Distributions for Bflav and Bcp candidates before 
vertex requirements: a) niES for -Bflav states; b) rriES for 
B" J/^ K°s,ip{2S)Kl,XciKl final states; and c) AS for 
the final state — > J/ip KI. In (a) and (b), the backgrounds 
are dominantly combinatorial. In (c) there are backgrounds 
from events containing a true J/^ but with a spurious K%. 
Other background comes from events in which no true J/i/> is 
present. 



a i?° or B^ when it decayed. This determination can- 
not be done perfectly. If the probability of an incorrect 
assignment is it;, an asymmetry that depends on the dif- 
ference between B^ and B^ tags will be reduced by a 
factor D = 1 — 2w, called the dilution. A neural network 
combining the outputs of algorithms that evaluate the 
characteristics of each event is used to take into account 
the correlations between the different sources of flavor 
information and to estimate and B^ mistag probabil- 
ities for each event. Based on these values and the source 
of flavor information, the event is tagged and assigned to 
one of five mutually exclusive tagging categories. The di- 
lution for each category is determined from the data, as 
described in Sec. VII. Grouping tags into categories, each 
with a relatively narrow range in mistag probability, in- 
creases the overall power of the tagging while simplifying 
the studies of systematic uncertainties. 

Events with an identified primary electron or muon 
and a kaon with the same charge, if present, are assigned 
to the Lepton category. Events with both an identified 
kaon and a low-momentum (soft) pion candidates with 
opposite charge and similar flight direction are assigned 
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to the KaonI category. Soft pion candidates from D*"*" 
decays are selected on the basis of their momentum and 
direction with respect to the thrust axis of Stag- Events 
with only an identified kaon are assigned to the KaonI 
or KaonI I category depending on the estimated mistag 
probability. Events with only a soft-pion candidate axe 
assigned to the KaonI I category as well. The remain- 
ing events are assigned to either the Inclusive or the 
UnTagged category based on the estimated mistag prob- 
ability. The UnTagged tagging category has a mistag rate 
set to 50%, and therefore does not provide tagging infor- 
mation. It does, however, increase the sensitivity to the 
decay-rate difference AF and allows the determination 
from the data of the detector charge asymmetries, as de- 
scribed in Sec. VII. This tagging algorithm is identical 
to that used in Ref. [5]. 

Wc consider separate mistag probabilities for B'^ and 
B° tags, w%f and w^, in each tagging category a. 
From these, we define the average mistag probability 
w" = +w'^^)/2 and the asymmetry in the mistag 
rates Aw" = — w'^ . A correlation between the av- 
erage mistag rate and the At uncertainty a/^t estimated 
evcnt-by-event (discussed in Sec. VI) is observed for kaon- 
based tags [2, 29]. For a At uncertainty less than 1.4 ps, 
this correlation is found to be approximately linear: 

w;«=<+<i„p,aAt . (29) 

All signal mistag parameters, Wq, w^^^^^, and Aw", are 
free in the global fit (11 in total since Wgi^'^p™ is assumed 
to be zero), and their results can be found in Table VIII 
in Sec. VIII. 



VI. DECAY-TIME MEASUREMENT AND At 
RESOLUTION FUNCTION 

The time interval At = tree ~ ^tag between the two 
B decays is calculated from the measured separation 
Az between the decay vertices of the reconstructed B^ec 
meson and the Stag meson along the ^;-axis, using the 
known boost of the T(45) resonance in the laboratory, 
/37 « 0.55, the beam-spot size, and the momentum of the 
fully reconstructed B meson. The method is the same as 
described in Sec. V in Ref. [29]. 

An estimated error (ta« on At is calculated for each 
event. This error accounts for uncertainties in the track 
parameters from the SVT and DCH hit resolution and 
from multiple scattering, for the beam-spot size, and for 
effects from the B-flight length transverse to the beam 
axis. However, it does not account for errors due to 
mistakes of the pattern recognition system, wrong as- 
sociations of tracks to vertices, misalignment within and 
between the tracking devices, inaccuracies in the model- 
ing of the amount of material in the tracking detectors, 
limitations in our knowledge of the beam-spot position, 
or uncertainty in the absolute z scale. Most of the ef- 



fects that are not explicitly accounted for in cta* sltc ab- 
sorbed in the At resolution function, described below. 
Remaining systematic uncertainties are discussed in de- 
tail in Sec. X. 

We use only those events in which the vertices of 
the -Brec and -Btag are successfully reconstructed and for 
which \At\ < 20ps and cta/ <1.4ps. The fraction of 
events in data satisfying these requirements is about 85%. 
From Monte Carlo simulation we find that the recon- 
struction efficiency does not depend on the true value 
of At. The r.m.s. Az resolution for 99.7% of the events 
used is about 160 /im (1.0 ps), and is dominated by the 
resolution of the -Btag vertex. 

To model the At resolution we use the sum of three 
Gaussian distributions (called core, tail and outlier com- 
ponents) with different means and widths: 

ftailhciSt] i^tailfAt, -StailfAt) + 

fouthG{St;6out,(^out) (30) 

where 

hG{St;S,a) = ^_e-(^*-^)'/(2a^) . (31) 
V2Tra 

Here St = At — Attruo represents the reconstruction er- 
ror and /core = 1 — /tail " /out- We incorporate the last 
Gaussian distribution in Eq. (30) without reference to 
(TAt since the outlier component is not expected to be 
well described by the estimated uncertainty. The widths 
of the first two Gaussian components are given by aAt 
multiplied by two independent scale factors, Score and 
'S'taii) to accommodate an overall underestimate {S > 1) 
or overestimate {S < 1) of the errors. The core and 
tail Gaussian distributions are allowed to have non-zero 
means ((5corofAt and i^taiiCAt) to account for residual bi- 
ases due to daughters of long-lived charm particles in- 
cluded in the Stag vertex. Separate means are used for 
the core distribution of each tagging category. These 
means are scaled by cta* to account for a correlation be- 
tween the mean of the St distribution and aAt [2, 29]. 
This correlation is found to be approximately linear for 
(TAt less than 1.4 ps. The non-zero means of the reso- 
lution function introduce an asymmetry into the other- 
wise symmetric At distributions. All other parameters 
of the resolution function are taken to be independent of 
the tagging category. We find that the three parameters 
describing the outlier Gaussian component are strongly 
correlated among themselves and with other resolution 
function parameters. Therefore, we fix the outlier bias 
Sout and width (Tout to Ops and 8ps, respectively, and 
vary them through a wide range to evaluate systematic 
uncertainties. The outlier Gaussian distribution accounts 
for less than 0.3% of the reconstructed vertices. 

In simulated events, we find no significant differences 
between the At resolution function of the Bflav, Bqpx°i 
and Bqpj^o samples. This is expected since the Stag ver- 
tex precision dominates the At resolution. Hence, the 
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same resolution function is used for all modes. Possible 
residual differences are taken into account in the evalua- 
tion of systematic errors described in Sec. X. 

The resulting signal resolution function is described 
by a total of 12 parameters, Scnr^, Sll^e"", S^^°f, .5^,^°°" 



^core ' ^core 



•^core"^^^^! '^core^^ , /tailj <5taili 'S'tail; /out) '^out) ""out) ten 

of which are free in the final fit. 

As a cross-check, we use an alternative resolution func- 
tion that is the sum of a single Gaussian distribution 
(centered at zero), the same Gaussian convolved with a 
one-sided exponential to describe the core and tail parts 
of the resolution function, and a single Gaussian distri- 
bution to describe the outlier component [2]. The expo- 
nential component is used to accommodate the bias due 
to tracks from charm decays originating from the -Btag- 
The exponential constant is scaled by aAt to account for 
the previously described correlation between the mean of 
the St distribution and aAt- In this case, each tagging 
category has a different core component fraction and ex- 
ponential constant. 



a B°, or a CP eigenstate. Treating K° and 
as if they were CP eigenstates introduces effects 
that are negligible on the scale of the statistical 
and systematic uncertainties of this analysis; 

iv) the decay-time measurement At and its estimated 
error cta*; 

v) a variable C used to assign the probability that the 
event is signal or background. Either ^ is rriES (for 
flavor eigenstates and CP eigenstates with K^) or 
it is AE (for CP eigenstates with if^). 

The likelihood function is built from time distributions 
that depend on whether the event is signal or any of a va- 
riety of backgrounds (together specified by the index j), 
on the tag category, on the tag flavor, and on the type 
of reconstructed final state. The contribution of a single 
event to the log-likelihood is 



log 



E 



•^"ec'(C)Ka^,rec(At,aAt) 



(32) 



VII. LIKELIHOOD FIT METHOD 

We perform a single, unbinned maximum-likelihood fit 

to all -Bfiav, BopjfQ , and Bfjpj^a samples. Each event is 
characterized by the following quantities: 

i) assigned tag category a G {Lepton, KaonI , 

Kaonll, Inclusive, UnTagged}; 
a) tag-flavor type "tag" € {Bt, Bt}, i.e., the tagging 
state is ostensibly a or B^, unless it is untagged; 
in) reconstructed event type "rec" G {Bf, Bf, CPK^, 
CPK^ }, i.e., the reconstructed state is ostensibly 



For a given reconstructed event type "rec" and tagging 
category a, J^^g'^iO gives the probability that the event 
belongs to the signal or any of the various backgrounds 
denoted by j. Each such component has its own prob- 
ability density function (PDF) W"^^ ^.^..(At, cta*), which 
depends as well on the particular tag flavor "tag" . This 
distribution is the convolution of a tagging-category- 
dependent time distribution H"^^ j.^^{AttTue) with a At 
resolution function TZ""'^ {St,<TAt) of the form given in 
Eq. (30), but with parameters that depend on the tag- 
ging category a and on the signal/background nature of 
the event j: 



/ + OC 
d(Attrue)7^"'^' {At - Attrue, CTAt )iftag',rec(Attrue) , (33) 
-oo 

where 

Fraiec(Attrue) = P^.ec {^tt'^' (1 " <a^)<,rec(Attrue) + T^^ftoi^/l^^^^ . (34) 



Here /itag,rec(^) represents the time dependence dN/dt 
given in Eqs. (17) - (21), with t = Attme- We indi- 
cate by the mistag fractions for category a and 

component j. The index "tag" denotes the opposite 
flavor to that given by "tag". For eveiits falling into 
tagging category UnTagged we define w^'^^— to be 1/2. 

The efficiency r^^g is the probability that an event whose 
signal/background nature is j and whose true tag fla- 
vor is "tag" will be assigned to category a, regardless 
of whether the flavor assigned is correct or not. The 
efficiency p^^^ is the probability that an event whose sig- 



nal/background nature is indicated by j and whose true 
reconstructed character is "rec" will, in fact, be recon- 
structed. For non-BB background sources, where the 
meaning of true "tag" and "rec" is ambiguous, this pro- 
vides an empirical description of the efficiencies as well 
as the mistag fractions. 
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A. PDF Normalization 

Every reconstructed event, whether signal or back- 
ground occurs at some time Atime^ so 



where k indexes the various conibinatorial (fc = prompt, 
non-prompt) background components, and 



a,k 
rec 



(40) 



/-l-oo 
d(Att,ue)/iLg,rec(Attrue) = 1 , (35) 
-oo 

for each value of "rec", "tag" and j. Moreover, ev- 
ery event is assigned to some tagging category (possibly 
UnTagged); thus 



(36) 



for each vahic of "tag" and j. It follows then that the 
normalization of H^g^^ .^^^{Attiue) is 

/ + 00 
d(Attrue)/?t"airec(Attrue) = 2p^,, . (37) 

« Lag 

In this analysis the nominal normalization of 

Wrag,roc(Ai>CrAt) iS the Same as i7t"g,rcc(A^truc), 

but fits with normalization in the interval [—20, 20] ps 
have been also performed as a cross-check to evaluate 
possible systematic effects. 



B. Signal and Background Chciracterization 

The function ^^'HC) in Eq. (32) describes the signal 
or background probability of observing a particular value 

of It satisfies 



/■Cmax , 



Oi,3 
rec 



(0 = 1 



(38) 



where [Cmin, Cmax] IS the range of niEs or AE values used 
for analysis. 

For Bflav and B^pj^o events, the tties shape is de- 
scribed with a single (faussian distribution for the sig- 
nal and an ARGUS parameterization for the back- 
ground [29] . Based on these fits, an event-by-event signal 
probability p"g(.(mEs) can be calculated for each tagging 
category a and sample "rec" . Since we do not expect sig- 
nal probability differences between and B^, the rriES 
fits are performed to Bf and Bf events together. The 
fits to B° i>{2S)K° and B" ^ xaK^ are performed 
without subdividing by tagging category, due to the lack 
of statistics and the high purity of the samples. We dis- 
tinguish three different background components: peaking 
background events, which have the same tties behavior as 
the signal; a zero-lifetime (prompt) combinatorial com- 
ponent; and a non-zero-lifetime (non-prompt) combina- 
torial background. The component fractions .F^'/(mEs) 
are then (j = sig, peak,fc) 



^ rec 



'(toes) 
^,:;,P^'^^(mEs) 



= [1 - /.".r^"^] p?,,(mEs), 



Jrcc Prcc(TOEs), 



•^rec ("^ES) = [l-p"ec("^Es)]/, 



a,k 
rec ' 



(39) 



The fraction f^^f'^'^^ of the signal Gaussian distribution 
is due to backgrounds that peak in the same regions as 
the signal, and is determined from Monte Carlo simula- 
tion [29]. The estimated contributions are (1.5 ± 0.6)%, 
(0.28±0.11)%, (1.8±0.6)%, (l+i^)%, and (3. 5± 1.4)% for 
the Bfiav, J/i^K", (if," ^ TT+TT-), J/ib K° {K'^ ^ 7r°7rO), 
^p{2S)Kg, and XaK^ channels, respectively. A com- 
mon peaking background fraction is assumed for all tag- 
ging categories within each decay mode. We also as- 
sume a common prompt fraction for all tagging cate- 
gories for each B(jpj^o decay channel. Since the -Bfiav 
sample is large and there are significant differences in the 
background levels for each tagging category, j"«j,P™™P* = 
j.a,prompt allowcd to depend on the tagging category. 

Note that the parameters of the T^^'^^^ (niEs) functions, 
determined from a set of separate unbinned maximum- 
likelihood fits to the mES distributions, are fixed in the 
global fit. 

For B(jpj^o events the background level is higher than 
it is for Bijpj^o , with significant non-combinatorial com- 
ponents [29] . A binned likelihood fit to the AE spectrum 
in the data is used to determine the relative amounts of 
signal and background from B — !■ J/tpX (e.g., J/ipK*) 
events and from events with misrcconstructed J/'ip 
£~^£~ candidates (non-J/'0 background). In these fits, 
the signal and B J/ipX background distributions 
are obtained from inclusive- J/?/) Monte Carlo samples, 
while the non- J/ip distribution is obtained from the J/^p 
dilepton-mass sideband. The Monte Carlo simulation is 
also used to evaluate the channels that contribute to the 
B J/tpX background. The fit is performed separately 
for if^ candidates reconstructed in the EMC and in the 
IFR, and for J/tp candidates reconstructed in the e+e~ 
and IJL~ modes, since there are differences in purity 
and background composition. Candidates reconstructed 
in both IFR and EMC are considered as belonging to the 
IFR category because of its better signal purity. The 
different inclusive- J/V^ backgrounds from Monte Carlo 
are then normalized to the J/V' background fraction ex- 
tracted from the AE fit in the data. The normaliza- 
tion to the data is performed separately for lepton-tagged 
and non-lepton-tagged events to account for the observed 
differences in flavor-tagging efficiencies between the J/ip 
sideband events and the Bflav and inclusive- J/^/; Monte 
Carlo events. In addition, some of the decay modes in 
the inclusive- J/'f/' background have CP content. The 
same PDFs are used to describe the AE shape for J/ij} 
candidates in the and e^e~ channels. However, 

different PDFs are used for K^s observed in the IFR and 
in the EMC. Separate AE PDFs are used for J/ipKl (sig- 
nal), J/tpKg background, J/tpX background (excluding 
J/tpK'^), and non- J/ijj background. 
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C. Efficiency Asymmetries 

For each signal or background j, the average recon- 
struction efficiencies p' = (p^^ + p'-^)/2, P^p^^o, and 

PcPK° absorbed into the fractions of reconstructed 
events falling into the different signal and background 
classes. In contrast, because all events fall into some tag- 
ging category (including UnTagged) , the average tagging 
efficiencies t"'^ = (j'^^ -\- 'r'^)/'^ are meaningful, and 
the fraction of untagged signal events plays an impor- 
tant role. The asymmetries in the efficiencies, 



7 7 



Bf 



'Bf 

Oi.j 



Bf 
Bt 



a.i 



+ T 



Bt 



(41) 



need to be determined precisely, because they might oth- 
erwise mimic fundamental asymmetries wc seek to mea- 
sure. In Appendix A we illustrate how the use of the un- 
tagged sample makes it possible to determine the asym- 
metries in the efficiencies. Note that asymmetries due 
to differences in the magnitudes of the decay amplitudes, 
I As/ 1 7^ \ and \ABt \ 7^ l^stl) cannot be distinguished 
from asymmetries in the efficiencies, and thus are ab- 
sorbed in the v and fi parameters. 

Wc determine the average tagging efficiencies t"^-' by 
counting the number of events falling into different tag- 
ging categories, without distinguishing where an event is 
signal or background (i.e., r"'^ =''■"), since for each tag- 
ging category a the j component dependence is absorbed 
into the fractions of events falling into the different sig- 
nal and background components. For signal events, the 
parameters v^^^ and are included as free parame- 

ters in the global fit, and are assumed to be the same 
for all peaking background sources. For i?+ peaking 
background components, v'p^^^ and ^"-p'*'^'^ are fixed to 
the values extracted from a previous unbinned maximum- 
likelihood fit to the tagged and untagged At distributions 
of data used as control samples, described in Sec. IV. 
For combinatorial background sources the v and n pa- 
rameters are neglected. 



D. Mistags and At Resolution Function 

For signal events, a common set of mistag and At reso- 
lution function parameters, independent of the particular 
fully reconstructed state, is assumed. This assumption is 
supported by Monte Carlo studies. 

Peaking backgrounds originating from B^ decays are 
assumed to have the same resolution function and mistag 
parameters as the signal. For i?"*" peaking backgrounds 
we assume the same resolution function as for signal, but 
the mistag parameters are fixed to the values extracted 
from the same maximum-likelihood fit to the B^ data 



used to extract the parameters yP'^'^^ and ^"■p'^'*'', as de- 
scribed above. 

For combinatorial background components (prompt 

and non-prompt components in the -Bflav and Bqpj^o 
samples and the non-J/ip background in the Bqpi^o sam- 
ple) we use an empirical description of the mistag proba- 
bilities and At resolution, allowing various intrinsic time 
dependences. The parameters Aw" and w^^p^ are fixed 
to zero, and the resolution model uses core and outlier 
Gaussian distributions. The fractions of prompt and non- 
prompt components and the lifetime of the non-prompt 
component in the non-J/ip background are fixed to the 
values obtained from an external fit to the time distribu- 
tion of the J/ip dilepton-mass sideband. 



E. EVee Parameters for the Nominal Fit 

The aim of the fit is to obtain simultaneously 
sgn(RcAcy)Ar/r, (Re Acy/|Acy |) Rez, and Imz, 

assuming rcp = 1. The parameters ImAcp/|Acp| and 
Am are also free in the fit to account for possible cor- 
relations and to provide an additional cross-check of the 
measurements. The average B^ lifetime tb = l/F is 
fixed to the PDG value, 1.542 ps [31]. As a cross-check 
we also perform fits allowing rcp and T to vary. All these 
physics parameters are, by construction, common to all 
samples and tagging categories, although the statistical 
power for determining each parameter comes from a par- 
ticular combination of samples or At dependences, as 
discussed in Sec. 11. 

The terms proportional to the real parts of the DCS- 
decay parameters are small since Kc Xsf and RcAg^ oc- 
cur only multiplied by other small parameters (see Ta- 
bles 11- IV), and are therefore neglected in the nomi- 
nal fit model. Fixing = 0.02, our best es- 
timate from \V*^Vj^^/V*ijV^^\, we fit for the parameter 
Im As/ /I As/ 1, and vary separately ImAgj/|Ag/|, keeping 

I^B/I = I^b/IIp/'jP- We do not require I Im As/ /As/ I < 1. 
Thus, there are two free parameters associated to DCS 
decays, plus one fixed magnitude. 

Wc treat Xst and Agj similarly. Since there is no inter- 
ference betweeen B° and B° semileptonic decays, we set 
^Bt = 0, Agj = for the Lepton tagging category. For 
the other tagging categories we assume common values 
of the DCS-dccay parameters. We assign a systematic 
error by varying \ABf/ABf\ and \ABt/ABt\ by 100% and 
scanning all possible combinations of the phases (Sec. X). 
With a larger data sample, direct determination of the 
DCS-dccay parameters might be advantageous. With the 
current sample, absorbing some of the variation into the 
systematic uncertainty suffices to prevent effects induced 
by DCS decays being misinterpreted as symmetry viola- 
tions. 

The total number of parameters that are free in the fit 

is 58, of which 36 parameterize the signal: physics pa- 
rameters (4), cross-check physics parameters (2), single 
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effective imaginary parts of the DCS-decay phases (4), 
resolution function (10), mistag probabilities (11), and 
differences in the fraction of 5° and -8° mesons that are 
tagged and reconstructed (5). The remaining 22 param- 
eters are used to model the combinatorial backgrounds: 
resolution function (3), mistag fractions (8), fractions of 
prompt components (9) and the effective lifetime of the 
non-prompt contributions (2). 

The At distributions, the asymmetries, the physics pa- 
rameters sgn(ReAcp)Ar/r, \q/p\, (Re Acp/|Acp|) Rez, 
and Imz and the cross-check parameter ImAcp/|Acp| 
were kept hidden until the analysis was finished. How- 
ever, the parameter Am, the residual At distributions 
and asymmetries, the statistical errors, and changes in 
the physics parameters due to changes in the analysis 
were not hidden. 



VIII. ANALYSIS RESULTS 

We extract the parameters sgn(Re Acp)Ar/r, \q/p\, 
(Re Acp/|Acp|) Rez, Imz, ImAcp/|Acp|, Am, the pa- 
rameters for DCS decays, the signal mistag probabili- 
ties, resolution-function and u and /i" parameters, and 
the empirical background parameters with the likelihood 
function described in Sec. VII. In Table VII we list the 
signal yields in each tagging category after vertex require- 
ments. The purities (estimated from the mEs fits for 
noii-Bopjft)^ samples and in the region |Ai?| < 10 MeV 
for Bqpko events) , averaged over tagging categories, are 
82%, 94%, and 55%, for B^av, ^cPifg , and Bf^pj^o candi- 
dates, respectively. The fitted signal mistag probabilities 
and resolution-function parameters are shown in Tables 
VIII and IX. The values of the asymmetries in recon- 
struction and tagging efficiencies are summarized in Ta- 
ble X. There is good agreement with the asymmetries 
extracted with the counting-based approach outlined in 
Appendix A. 

TABLE VII: Signal event yields after vertex requirements, 
obtained from the mES fits for the Bflav and Bqpi^o samples. 
For the B/jpj^o sample, the signal yields are obtained using 
the signal fractions determined from the fit to the AE dis- 
tributions, and are quoted for events satisfying \AE\ < 10 
MeV. 









Bcpro 


Tag 


B° Tot 


B" Tot 


B° W Tot 


Lepton 
Kaon I 
Kaon II 
Inclusive 
UnTagged 


1478 1419 2897 
2665 2672 5337 
3183 2976 6159 
3197 3014 6211 
10423 


96 98 194 
154 175 329 

181 188 369 
184 172 356 
585 


35 35 70 
74 65 139 

85 66 151 
78 72 150 
260 



The values of the parameters sgn(Rc Acp)Ar/r, \q/p\, 
(Re Acp/|Acp|) Rez, and Imz extracted from the fits are 



TABLE VIII: Average tagging efficiencies after vertex require- 
ments and signal mistag parameters for each tagging category 
a as extracted from the maximum-likelihood fit that allows 
for CPT violation. Uncertainties are statistical only. 



Tagging 


T«(%) 


<■"«(%) 


a,sig 
^slope 


Aw«'='s(%) 


category 










Lepton 


9.4 ±0.2 


2.6 ±0.7 


(fixed) 


-1.2 ± 1.2 


Kaon I 


17.2 ±0.3 


2.0 ± 2.0 


0.13 ±0.04 


-2.7 ±1.3 


Kaon II 


19.9 ± 0.3 


15.9 ±2.4 


0.07 ±0.04 


-4.2 ±1.3 


Inclusive 


19.9 ± 0.3 


26.5 ± 2.5 


0.07 ±0.04 


-2.9 ±1.3 


UnTagged 


33.6 ± 0.6 


50 (fixed) 


(fixed) 


(fixed) 



TABLE IX: Signal At resolution function parameters as ex- 
tracted from the maximum-likelihood fit that allows for CPT 
violation. Uncertainties are statistical only. 



Parameter 


Fitted value 


Parameter 


Fitted value 


'5'corc 


1.25 ± 0.04 


"Stail 


5.7 ±0.8 


eLepton 

Ocore 


0.02 ±0.07 


5tail 


-1.5 ±0.5 


cKaonl 
^core 


-0.27 ±0.05 


/tail 


0.034 ±0.010 


cKaonll 

'^core 


-0.32 ± 0.04 


Cout 


8 ps (fixed) 


cinclusive 
^corc 


-0.30 ± 0.04 


^out 


ps (fixed) 


cUnTagged 
Ocore 


-0.28 ± 0.03 


/out 


0.0003 ± 0.0012 



TABLE X: Values of the signal differences m recon- 

struction {v''^^) and tagging (/i™'^'^) efficiencies as extracted 
from the maximum-likelihood fit that allows for CPT viola- 
tion. The results are compared with those obtained with a 
counting-based method described in Appendix A. 



Parameter 


Nominal fit 


Counting-based method 


^Lepton, sig 
Kaonl.sig 
KaonII,sig 

A* 

Inclusive, sig 

A* 


0.011 ± 0.008 
0.024 ± 0.022 

-0.022 ±0.017 
0.014 ±0.016 
0.014 ± 0.016 


0.007 ± 0.008 
0.029 ± 0.042 

-0.022 ± 0.029 
0.004 ± 0.027 
0.025 ± 0.027 



given in Table XI. The fitted ImAcp/|Acp|, Am, and 
effective DCS-decay parameters are also indicated. All 
these results can be compared to those obtained when 
the fit is repeated assuming CPT invariance. The change 
in the effective DCS-decay parameters between the two 
fits is due to the large correlation of these parameters 
with the CPT- violating parameter Im z. The fitted value 
of Am agrees with recent B-factory measurements [1- 
4], and remains unchanged between the two fits. The 
fit result for Im Acp / 1 A^p | when we assume CPT invari- 
ance agrees with our sin 2/3 measurement based on the 
same data set [5]. When we allow for CPT violation, 
ImAcp/|Acp| increases by ±0.011, equal to 15% of the 
statistical uncertainty on ImAcp/|Acp|, which is con- 
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sistcnt with the statistical correlations observed in the 
fit with z free. The correlation coefficients among all 
physics and cross-check physics parameters arc shown in 
Table XII. The largest observed correlation (17%) ap- 
pears between Imz and ImXcp /\\cp\- Table XIII shows 
the largest statistical correlations of the physics param- 
eters with any other free parameter in the fit. Note that 
the variables \q/p\ and i/^'^ are significantly correlated, as 
arc Im z and the DCS-decay parameters. We do not eval- 
uate the full systematic errors for Am and Im Acp / 1 \cp \ 
so these measurements do not supersede previous BABMi 
measurements for these quantities. 

TABLE XI: Physics parameters extracted from the 
maximum-likelihood fits both allowing for CPT violation and 
excluding it. The free DCS-decay parameters are also indi- 
cated. Errors are statistical only. 



Parameter 


Fit with z free 


Fit with z = 


sgn(ReAcp)Ar/r 


-0.008 ±0.037 


-0.009 ± 0.037 


\1/P\ 


1.029 ±0.013 


1.029 ± 0.013 


(ReAcp/|Acp|)Rez 


0.014 ±0.035 




Imz 


0.038 ± 0.029 




Arn ( ps^"') 


0.521 ± 0.008 


0.521 ± 0.008 


Ini A(^/|Acx• 


0.752 = ().0()7 


().7ii ^ 0.1)67 


Im Ast/|Ast 


1.5 ± 1.2 


0.5 ±1.0 


Im Ai3t/|Asi 


-0.1 ± 1.2 


0.8 ±1.0 


Im As// Ac/' 


2.3 ± 1.1 


1.4 ±0.9 


Im Ac// Ac/i 


-O.G ^ 1.1 


O.i ± 0.9 



TABLE XIII: The largest correlations of each physics param- 
eter with other free parameters of the maximum-likelihood 
fit. 



Physics parameter 


Parameter 


Correlation (%) 


Am 


Lepton.sig 


-20.1 




/tail 


18.7 




5'tail 


-15.4 


sgn(ReAcp)Ar/r 


\<1/P\ 


11.0 






65.1 






-22.5 




Lepton,sig 

A* 


22.4 




^^KaonI,sig 


-22.4 




^^Inclusive.sig 


-15.5 




Kaonl.sig 


13.9 






-13.5 




sgn(ReAcp)Ar/r 


11.0 


ImAcy/IAcpl 


Imz 


17.4 




Im Ab*/ Ab* 


14.4 




Im As/ /I As/ 1 


13.6 




Rez 


-10.9 


(ReAcp/IAcyDRez 


ImXcp/\Xcp\ 


-10.9 


Imz 


Im Ast/ Ast 


61.6 




ImAs//|As/| 


57.7 




Im Ast/ Ast 


-56.6 




ImAsf/IAs/l 


-54.0 




Im Acf/ Acp| 


17.4 






11.0 



TABLE XII: Correlation (in %) among all the physics param- 
eters extracted from the simultaneous maximum-likelihood fit 
to the Bflav and Bcp samples. 



Parameter 


Parameter 


Correlation (%) 






z free 


z = 


Am 


sgn(ReAGp)Ar/r 


-1.3 


-0.9 




\1/P\ 


-2.8 


-2.8 




ImAcy/IAcpl 


-5.6 


-5.3 




(ReAcp/IAcyDRez 


7.0 






Imz 


-0.2 




sgn(ReAGp)Ar/r 


\1/P\ 


11.0 


10.8 




ImAcy/|Acy 1 


0.4 


0.2 




(ReAcp/IAcyDRez 


-7.9 






Imz 


-1.8 




Il/Pl 


Im Xcp/\\cp\ 


-1.0 


-1.5 




(ReAcp/|A(y|)Rez 


-2.4 






Imz 


-1.1 




ImAcy/IAcyl 


(ReAcp/IAcyDRez 


-10.9 






Imz 


17.4 




(ReAcp/IAcpDRez 


Imz 


-3.4 





Figs. 2 and 3 show the At distributions of events con- 
fined to the signal region, defined as thes > 5.27 GeV/c^ 
for the -Bflav and B^pj^a samples, and | A_E| < 10 McV for 
the Bqpj^o sample. The points correspond to data. The 
curves correspond to the projections of the likelihood fit 
allowing for CPT violation, weighted by the appropriate 
relative amounts of signal and background. The back- 
ground contribution is indicated by the shaded area. 



IX. CROSS-CHECKS AND VALIDATION 
STUDIES 

We use data and Monte Carlo samples to perform val- 
idation studies of the analysis technique. The Monte 
Carlo tests include studies with parameterized fast Monte 
Carlo as well as full GEANT4-simulated samples. Checks 
with data arc performed with control samples, where no 
AT and CP-, T-, and CPT- violating effects are expected. 
Other checks are made by analyzing the actual data sam- 
ple, but using alternative tagging, vertexing, and fitting 
configurations. 
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FIG. 2: The At distributions for (a) mixed and (b) unmixed 
Bflav events with a tag or with a tag in the signal 
region, tties > 5.27 GeV/c^. The solid (dashed) curves rep- 
resent the fit projection in At based on the individual sig- 
nal and background probabilities and the event-by-event At 
uncertainty for B° tags. The shaded area shows the 

background contribution to the distributions. 



A. Monte Carlo Simulation Studies 

A test of the fitting procedure is performed with pa- 
rameterized Monte Carlo simulations consisting of 300 
experiments generated vifith a sample size and compo- 
sition corresponding to that of the data. The mistag 
probabilities and At distributions are generated accord- 
ing to the model used in the likelihood function. The 
physics parameters are generated according to the values 
found in the data [34] . The nominal fit is then performed 
on each of these experiments. Each experiment uses the 
set of Toes (A-B) and aAt values observed in the non-i^2 
{K1) sample. The r.m.s. spread of the residual distri- 
butions for all physics parameters (where the residual is 
defined as the difference between the fitted and generated 
values) is found to be consistent, within 10%, with the 
mean (Gaussian) statistical errors reported by the fits. 
Moreover, it has been verified using these experiments 
that the asymmetric 68% and 90% confidence-level inter- 
vals obtained from the fits provide the correct statistical 
coverage. 

In all cases, the mean values of the residual distribu- 
tions are consistent with no measurement bias. A sys- 
tematic error due to the limited precision of this study is 
assigned to each physics parameter. The statistical errors 
on all the physics parameters (Table XI) and the calcu- 
lated correlation coefficients among them (Tables XII and 
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FIG. 3: The At distributions for (a) -BcPifg ^-nd (b) B^pj^o 
events with a tag or with a tag in the signal region, 
»nES > 5.27 G^/c^ for Bop^a candidates and |A_B| < 10 
MeV for -Bq^^o events. The solid (dashed) curves represent 
the fit projection in At based on the individual signal and 
background probabilities and the event-by-event At uncer- 
tainty for B° (B°) tags. The shaded area shows the back- 
ground contribution to the distributions. 



XIII), extracted from the fit are consistent with the range 
of values obtained from these experiments. We find that 
24% of the fits result in a value of the log-likelihood that 
is greater (better) than that found in data. 

In addition, samples of signal and background Monte 
Carlo events generated with a full detector simulation are 
used to validate the measurement. The largest samples 
are generated with AF/F, \q/p\ — 1, and z all equal to 
zero, but additional samples are also produced with rela- 
tively large values of these parameters. Other values (in- 
cluding those measured in the data) are generated with 
reweighting techniques. The signal Monte Carlo events 
are split into samples whose size and proportions of -Bflav, 
Bqpj^o , and Bqpj^o are similar to those of the actual 
data set. To check whether the selection criteria or the 
analysis and fitting procedures introduce any bias in the 
measurements, the fit (to signal alone) is then carried out 
on these experiments, allowing for CPT violation. The 
small combinatorial background in these signal samples is 
suppressed by restricting the fit to the events in the signal 
region. Fits to a sample without background, using the 
true At distribution and true tagging information, are 
also performed. The means of the residual distributions 
from all these experiments for all the physics parameters 
are consistent with zero, confirming that there is no mea- 
surement bias. The r.m.s. spreads are consistent with the 
average reported errors. A systematic error is assigned 
to each physics parameter corresponding to the limited 
Monte Carlo statistics for this test. 
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The effect of baekgrounds is evaluated by adding an 
appropriate fraction of background events to the signal 
Monte Carlo sample and performing the fit. The Bcp 
background samples are obtained either from simulated 
B J/tl>X events or AE sidebands in data, while the 
Bfiav backgrounds are obtained from generic BB Monte 
Carlo. We find no evidence for bias in any of the physics 
parameters. 



B. Cross-checks with Data 

We fit subsamples defined by tagging category or data 
taking period. Fits using only the B^ — > D^*^~X~^ 
or B° J/tp K*°{K+ Tr~) channels for Ba^v, and only 
BcPK% or only B^pj^o for Bcp are also performed. We 
find no statistically significant differences in the results 
for the different subsets. We also vary the maximum 
allowed values of \At\ between 5 and 30 ps, and of a At 
between 0.6 and 2.2 ps. Again, we do not find statistically 
significant changes in the physics parameters. 

In order to verify that the results are stable under vari- 
ation of the vertex algorithm used in the measurement 
of At, we use alternative (less powerful) methods, de- 
scribed in Sec. VIII.C.5 in Ref. [29]. To reduce statistical 
fluctuations due to different events being selected, the 
comparison between the alternative and nominal meth- 
ods is performed using only the events accepted by both 
methods. Observed variations are small compared with 
the systematic error assigned to the resolution function 
(see Sec. X). 

The stability of the results under variation of the tag- 
ging algorithm is studied by repeating the fit using the 
tagging algorithm described in Sec. IV in Ref. [29]. The 
algorithm used in that analysis has an effective tagging 
efficiency Q = X)a ^ 2w")^ about 7% lower than 
the one used here. The variations observed in the physics 
parameters are consistent with the statistical differences. 

The average B'^ lifetime is fixed in the nominal fit to 
the PDG value [31]. This value is obtained by averag- 
ing measurements based on fiavor-eigenstate samples and 
by assuming negligible effects from AF/F, \q/p\, and z. 
Measurements that do not use tagged events are largely 
insensitive to \q/p\ and z, but would be affected, at sec- 
ond order, by a non-zero value of AF/F, as discussed in 
Sec. II. Therefore we do not expect sizeable effects from 
the fixed average lifetime. However, to check the con- 
sistency of the result, the fit is repeated with the average 
lifetime left free. The resulting tb is about one stan- 
dard deviation below the nominal value assumed in our 
analysis, taking into account the statistical error from 
the fit and the present tb uncertainty. As described in 
Sec. X C, a systematic error is assigned using the varia- 
tion of each physics parameter when the fit is repeated 
with Tb fixed to the value obtained when it is floated, 
which corresponds to a change of about twice the present 
PDG error (±0.032 ps). 

Similarly, flts with rcp free have been performed. The 



resulting rcp value is consistent with unity (the fixed 
nominal value) within one standard deviation (statistical 
only). As described in Sec. X C, systematic errors due to 
fixing rcp at unity are set by changing rcp by twice the 
statistical uncertainty determined by leaving it free in the 
fit (±10%). The resulting variation in each parameter is 
taken as the systematic error. 

The robustness of the fit is also tested by modifying the 
nominal PDF normalization, as described by Eq. (37), 
so that the analysis is insensitive to the relative number 
of B^ and B^ tagged events. As a consequence, the sta- 
tistical error on \q/p\ is dramatically increased, since the 
sensitivity to this parameter comes largely from the dif- 
ferences in time-integrated and B^ rates. In addition, 
the fit is also performed assuming an independent set of 
resolution function parameters for each tagging category. 
In all cases the results are consistent with the nominal 
fit results. Finally, the tagging efficiencies r" are alter- 
natively determined for each sample (Bflav, B^pj^o and 
BcPK° ) separately, rather than using a common estimate 
from tlie .Bflav sample, as in the nominal fit. The changes 
in the values of the physics parameters are negligible. 

Control samples in data from B~^ decays (treated in 
a way analogous to that described in Sec. IV) are also 
used to validate the analysis technique, since in these 
samples wc expect zero values for AF/F, \q/p\ — 1 and 
z. For the Bflav sample we use the B+ —f D"tt+ ,D*°Tr~^ 
decay channels, and for the Bcp sample the decays of 
charged- i? mesons to charmonium plus a charged K or 
K* (see Table VI). The check is performed by fixing 
Am = and \q/p\ = 1 in the Bflav sample, and assum- 
ing maximal mixing (Am = 0.489 ps~^ [31]) in the Bcp 
sample, and fitting for Im Acp / 1 Acp | , sgn(Re Acp)Ar/F, 
(Re Acp/|Acp|) Rcz and Imz. No statistically significant 
deviations from zero are observed. 



X. SYSTEMATIC UNCERTAINTIES 

We estimate systematic uncertainties with studies per- 
formed on both data and Monte Carlo simulation sam- 
ples. A summary of the sources of non-negligible uncer- 
tainties is shown in Table XIV. In the following, the in- 
dividual contributions are referenced by the lettered lines 
in the Table. 



A. LikeHhood Fit Procedure 

Several sources of systematic uncertainties due to the 
likelihood fit procedure are considered. We include the 
results from the tests performed using the parameter- 
ized Monte Carlo sample (a) and the full GEANT4 signal 
Monte Carlo sample (b), as described in Sec. IX A. No 
statistically significant bias (mean of the residual distri- 
butions) is observed. Thus, we assign a systematic error 
equal to the statistical uncertainty on the bias. No cor- 
rections are applied to the central values extracted from 
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TABLE XIV: Summary of systematic uncertainties on the measurements of sgn(Re Acp)Ar/r, \q/p\i (Re Acy/jAcpl) Rez, and 
Imz. 



Systematics source 


sgn(RcAG7')Ar/r 


\q/p\ 


(ReAGp/|Acp|)Rez 


Imz 


Likelihood fit procedure 


(a) Parameterized MC test 


0.003 


0.001 


0.003 


0.003 


(b) GEANT4-simulation test 


0.005 


0.007 


0.004 


0.016 


At resolution function 


(c) Resolution function parameterization 


0.007 


0.001 


0.008 


0.003 


(d) z scale and boost 


0.003 


0.001 


0.002 


< 0.001 


(e) Beam spot 


0.008 


0.002 


0.001 


0.011 


(f) SVT alignment 


0.006 


0.001 


0.001 


0.011 


(g) Outliers 


0.002 


< 0.001 


< 0.001 


< 0.001 


Signal properties 


(h) Average B° lifetime 


0.004 


0.001 


0.004 


< 0.001 


(i) Direct CP violation 


0.002 


0.004 


0.001 


0.003 


(j) DCS decays 


0.008 


0.004 


0.032 


0.006 


(k) Residual charge asymmetries 


0.005 


0.006 


0.004 


0.006 


Background properties 


(1) Signal probability 


0.002 


0.001 


0.002 


0.001 


(m) Fraction of peaking background 


< 0.001 


< 0.001 


0.004 


< 0.001 


(n) At structure 


0.002 


0.001 


0.001 


0.001 


(o) AF/CP/T/CPr/Mixing/DCS content 


0.001 


0.002 


0.002 


< 0.001 


(p) Residual charge asymmetry 


< 0.001 


0.001 


< 0.001 


< 0.001 


(q) /sT^-specific systematic errors 


0.004 


< 0.001 


0.004 


0.003 


Total systematic uncertainties 


0.018 


0.011 


0.034 


0.025 



the fit to the data. Note that the GEANT4 contribu- 
tion accounts for residual differences between the Sfiav, 
Bqpxo, and B^p^o samples in the mistag probability, 
resolution function, and z/ and /i" parameters. It also 
includes residual differences in At resolution for correct 
and wrong tags. 

We also consider the impact on the measured physics 
parameters of normalizing the time-dependent PDFs to 
the full interval — oo < < oo. The effect is evaluated 
by repeating the fit using a normalization in the range 
defined by the At cut. Finally, the fixed tagging efficien- 
cies are varied within their statistical uncertainties. The 
two contributions are negligible. 



B. At Resolution Function 

The At resolution model used in the analysis, con- 
sisting of the sum of three Gaussian distributions, is 
expected to be flexible enough to represent the exper- 
imental resolutions. To assign a systematic error for 
this assumption we use the alternative model described 
in Sec. VI, with a Gaussian distribution plus the same 
Gaussian convolved with one exponential function, for 
both signal and background. The results for all physics 



parameters obtained from the two resolution models arc 
consistent and we assign the difference of central values 
as a systematic uncertainty (c). 

In addition, a number of parameters that are inher- 
ent to the determination of At are varied according to 
known uncertainties. The PEP-II boost, estimated from 
the beam energies, has an uncertainty of 0.1% [27]. The 
absolute 2;-scale uncertainty is evaluated to be less than 
0.4%. This estimate is obtained by measuring the beam 
pipe dimensions with scattered protons and comparing 
to optical survey data. Therefore, the boost and z-scale 
systematic uncertainties are evaluated conservatively by 
varying by ±0.6% the reconstructed At and a At (d). As 
the beam spot is much smaller in the vertical than in 
the horizontal dimension, its vertical position and size is 
more relevant in the vertex fits. Hence the uncertainty on 
the position and size of the beam spot used in the vertex 
fits is taken into account by changing the vertical posi- 
tion by up to 40 /xm and increasing the vertical size from 
10 yum to 60 ^m (e) . Finally, the systematic uncertainty 
due to possible SVT internal misalignment is evaluated 
by applying a number of possible misalignment scenarios 
to a sample of simulated events and comparing the val- 
ues of the fitted physics parameters from these samples 
to the case of perfect alignment (f). 
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Fixing the width and bias of the outlier Gaussian dis- 
tribution in the resolution function to 8 ps and ps, re- 
spectively, is a potential source of bias. To estimate the 
corresponding systematic uncertainty we add in quadra- 
ture the variation observed in the physics parameters 
when the bias changes by ±5 ps, the width varies between 
6ps and 12 ps, and the outlier distribution is assumed to 
be flat (g). 



C. Signal Properties 

As described in Sec. IX B, the uncertainty from fixing 
the average i?° lifetime is evaluated by changing its cen- 
tral value by (±0.032 ps)-^ (h), twice the PDG error [31]. 
Possible direct CP violation in the Bqp sample is taken 
into account by varying rcp by ±10% (i). 

Systematic uncertainties related to DCS decays arise 
because we fix the real parts of A^t , A^^ , As/ , and A^ to 
zero. In order to evaluate this contribution, we generate 
samples of parameterized Monte Carlo samples tuned to 
the data sample, scanning the DCS-decay phases over 
their full allowed range (0-27r) and assuming a single 
hadronic decay channel contributing to the Stag and to 
the i?flav Samples arc generated with values of lAst/Astl 
and \ABf/ABf \ equal to and 0.04, corresponding to 
100% variation of the value 0.02 used in the nominal 
fit. For the Lepton tagging category, dominated by 
semileptonic B decays, wc assume Xbl to be zero. While 
the ratio of CKM matrix elements leads to the nominal 
value \ABt/ABt\ = |^b//Ab/| = 0.02, this is not a reli- 
able estimate for any single decay mode. Examination 
of the DCS charmed-meson decay I>° K'^tt~ shows 
good agreement with expectations from CKM matrix el- 
ements, albeit with large uncertainties, but the singly- 
CKM-suppressed decays D° tt+tt^ and K^K^ 
show deviations as large as a factor of two. However, 
when we sum over many channels, as we do here both for 
tagging states and for flavor eigenstates, quark-level pre- 
dictions arc much more reliable than they arc for a single 
channel. Allowing for 100% variation from the nominal 
value of 0.02 is thus conservative. 

Using the fit results from all these samples, we de- 
termine the offsets with respect to the generated value 
and its statistical uncertainty, for a complete sampling 
of DCS-decay phases. The systematic error assigned is 
the largest value among all configurations (j). This is 
the dominant source of systematic uncertainty for the 
measurement of (Re Acp/|Acp|) Rez and is due primarily 
to the influence of DCS decays in the tagging- _B meson. 
The effect of using a single effective channel for the fla- 
vor and all tagging category states has been estimated 
by splitting the Bflav and Stag samples generated with 
the parameterized Monte Carlo into equally sized sub- 
samples. For the different combinations of DCS-decay 
phases, the observed offset is about the average of the 
biases obtained using the single effective channel. There- 
fore, the largest offset among all conflgurations is smaller 



than that observed for a single channel. This shows that 
our prescription to describe the effects from DCS decays 
and to assign the systematic uncertainties assuming a 
single effective channel is conservative. 

Charge asymmetries induced by a difference in the de- 
tector response for positive and negative tracks are in- 
cluded in the PDF and extracted together with the other 
parameters from the time-dependent analysis. Thus, 
they do not contribute to the systematic error, but rather 
are incorporated into the statistical error at a level deter- 
mined by the size of the Bflav data sample. Nevertheless, 
in order to account for any possible residual effect, wc 
assign a systematic uncertainty as follows. We rerun the 
B reconstruction, vertex-flnding, and tagging algorithms 
after removing randomly and uniformly (no momentum 
or angular dependence) 5% of positive and, separately, 
negative tracks in the full Monte Carlo sample. This 
value of 5% is on average more than a factor of three 
larger than the precision with which the parameters v^^^ 
and have been measured in the data. Half the 

difference between the results obtained for positive and 
negative tracks is assigned as a systematic error (k). 



D. Background Properties 

The event- by- event signal probability p^(.(mEs) for 
Sfiav and Bqpj^o samples is fixed to the values obtained 
from the jdes Sts. We compare the results from the 
nominal fit to the values obtained by varying all the 
mES distribution parameters by ±1<t, taking into account 
their correlations. This is performed simultaneously for 
all tagging categories, and independently for the Sflav 
and Bopj^n samples. As an alternative, we also use a 
flat signal probability distribution: events belonging to 
the sideband region (toes<5.27 GeV/c^) are assigned a 
signal probability of zero, while we give a signal prob- 
ability equal to the purity of the corresponding sample 
to signal region events (mES>5.27 GeV/c^). The differ- 
ences among fitted physical parameters with respect to 
the default method are found to be consistent. We de- 
termine the systematic error due to this parameterization 
by varying the signal probability by its statistical error. 
The final systematic error is taken to be the larger of 
the one-sigma variations found for the two methods (1). 
The uncertainty on the fraction of peaking background 
is estimated by varying the fractions according to their 
uncertainties separately for the i?fiav sample and each 
Bqpj^o decay mode (m) . The effective rjcF of the Bqpj^o 
peaking background, assumed to be zero in the nominal 
fit, is also varied between +1 and —1 and the variations 
induced are negligible. 

Another source of systematic uncertainty is the as- 
sumption that the At behavior of the combinatorial back- 
ground in the toes sideband region is the same as it is in 
the signal region. However, the background composition 
changes gradually as a function of toeSj since the fraction 
due to continuum production slowly decreases as toes in- 
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creases. To study the impact of variable At behavior over 
the JTiEs range, we vary the lower edge of the mgg distri- 
butions from 5.20GeV/c^ to 5.27GeV/c^, simultaneously 
for the Bflav and Bqpj^o samples, observing good stabil- 
ity in the results. We also split the sideband region in 
seven equal slices each 10 MeV/c^ wide and repeat the 
fit in each of these slices. The results obtained for all 
physics parameters and m-ES slices are then linearly ex- 
trapolated to the i?-mass signal region. The quadratic 
sum of the extrapolation and the error on it is assigned 
as a systematic uncertainty (n). 

As described in Sec. VII, the likelihood fit assumes 
that there are no effects of AT, CP, T or CPT viola- 
tion, mixing, and DCS decays in the combinatorial back- 
ground components (-Bflav and Bqpj^o samples) and in 
the non-J/t/j background (iJ^y^o sample). To evaluate 
the effect of this assumption, we repeat the fit assum- 
ing for the background non-zero values of AF, \q/p\ — 1, 
z, ImAcp/|Acp|, and Am, and varying rjcp of the back- 
ground by ±1. The check is performed by introducing an 
independent set of physics parameters in the PDF and 
assuming maximal mixing and CP violation (Am and 
lm\cp/\Xcp\ fixed to 0.489ps-i [31] and 0.75 [5], re- 
spectively). DCS-decay effects are included by assuming 
the maximal values (0.04) of lAst/Astl and |As//Ab/|, 
and scanning all the possible values of the and 
phases for Sflav and Stag- The systematic uncertainty is 
evaluated simultaneously for all these sources (o). 

The systematic errors due to the _B+ decay rate 
are evaluated by varying its value by the PDG uncer- 
tainty [31]. The effect is negligible. The B+ mistags 
and the differences in the fraction of _B+ and B~ mesons 
that are tagged and reconstructed arc varied according 
to their statistical errors as obtained from the fit to the 
5+ data. These errors are found to be negligible. 

Uncertainties from charge asymmetries in combinato- 
rial background components (neglected in the nominal 
fit) are evaluated by repeating the fit with a new set of 
u and y^" parameters. The measured values of i/ and ^" 
are found to be compatible with zero and the variation 
of the physical parameters with respect to the nominal 
fit is assigned as a systematic error (p). 

For the Bfjpj^o channel, the signal and non- J/ip back- 
ground fractions are varied according to their statisti- 
cal uncertainties, obtained from the fit to the AE dis- 
tribution. We also vary background parameters, includ- 
ing the J/ipX branching fractions, the assumed rjcp, the 
AE shape, and the fraction and effective lifetime of the 
prompt and non-prompt non-J/^J components. The dif- 
ferences observed between data and Monte Carlo simula- 
tion for the Kl angular resolution and for the fractions of 
B° J/tpK^ events reconstructed in the EMC and IFR 
are used to evaluate a systematic uncertainty due to the 
simulation of the reconstruction. Finally, an addi- 
tional contribution is assigned to the correction applied 
to Lepton events due to the observed differences in flavor 
tagging efficiencies in the J/ip sideband relative to i?flav 
and inclusive J/tp Monte Carlo samples. Conservatively, 



this error is evaluated by comparing the fit results with 
and without the correction. The total B^pj^o -specific 
systematic error is evaluated by taking the quadratic sum 
of the individual contributions (q). 

E. Summary of Systematic Uncertainties 

All individual systematic contributions described 
above and summarized in Table XIV are added in quadra- 
ture. The dominant source of systematic error in the 
measurement of (Re Acp/|Acp|) Rcz is due to our lim- 
ited knowledge of the DCS decays, which also contributes 
significantly to the uncertainties on the other measure- 
ments. The limited Monte Carlo sample size is a dom- 
inant source of systematic error for \q/p\, Imz, and to 
a lesser extent for sgn(Re Acp)Ar/r. Residual charge 
asymmetries, mainly due to limited simulation statis- 
tics, dominate the systematic error on \q/p\. Our limited 
knowledge of the beam spot and SVT alignment reflects 
significantly on Imz and sgn(Rc Acp)Ar/r. The system- 
atic error on sgn(Re Acp)Ar/r receives a non- negligible 
contribution from our incomplete understanding of the 
resolution function. 

The systematic uncertainties on sgn(Re Acp) AF/F and 
\q/p\ when CPT invariancc is assumed are evaluated sim- 
ilarly, and found to be consistent, within the statistical 
fluctuations of the Monte Carlo simulation, with those 
found for the analysis when CPT violation is allowed. 



XI. SUMMARY AND DISCUSSION OF 
RESULTS 

The conventional analyses of mixing and CP violation 
in the neutral- _B-meson system neglect possible contribu- 
tions from several sources that are expected to be small. 
These include the difference of the decay rates of the 
two neutral-B-meson mass eigenstates, the CP- and T- 
violating quantity \q/p\ — l, and potential CPT violation. 
To measure or extract limits on these quantities requires 
the full expressions for time dependence in mixing and 
CP violation and consideration of systematic effects that 
might mimic the fundamental asymmetries we seek to 
measure. Such systematic effects could be induced by 
detector charge asymmetries, different rcsohition func- 
tions for positive and negative At, and DCS decays for 
both fully reconstructed final flavor states and nonlep- 
tonic tagging states. 

A hmit on the decay-rate difference of |AF/F| < 80% 
at 95% confidc;ncc;-level was obtained by CLEO [9] using 
the time-integrated mixing parameter Xd and the mass 
difference Am extracted under the assumption AF = 0. 
Using Z° decays, DELPHI [10] has recently performed 
a time-dependent study of semileptonic B decays inclu- 
sively reconstructed. Assuming no CP, T, or CPT vio- 
lation in mixing, they quote the limit |AF/F] < 18% at 
95% confidence-level. 
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Both \q/p\ and Imz were measured by OPAL [20], us- 
ing Z° decays to bb pairs and assuming AF = 0. Neutral- 
B-meson oscillations were studied by observing a single 
lepton indicative of a -B decay and the jet-charge asso- 
ciated with both the jet containing the lepton and the 
other jet. Because the multiparticle final states provide 
essentially uncorrelated B mesons, the issue of DCS de- 
cays is obviated. The results were Re cb = 0.006 ± 
0.010 ±0.006, equivalent to \q/p\ = 0.988 ±0.020 ±0.012, 
and Imes = -0.020 ± 0.016 ± 0.006, equivalent to 
Imz = 0.040 ± 0.032 ± 0.012. Combining the earlier 
\q/p\ measurements, all obtained assuming AF = 0, gives 
\q/p\ = 0.9993 ± 0.0064 [32]. Belle has used dilepton 
events to obtain limits on CPT violation [33]. Assuming 
that AF = and that CP violation in mixing can be 
ignored, they find Re cos 6* = - Re z = 0.00 ± 0.12 ± 0.02 
and Imcos6» = - Imz = 0.03 ± 0.01 ± 0.03. 

Our analysis of approximately 31,000 fully recon- 
structed flavor eigenstates and 2,600 CP eigenstates sets 
new limits on the difference of decay rates of mesons, 
and on the CP, T, and CPT violation intrinsic to B°B° 
mixing. The six independent parameters governing oscil- 
lations (Am, AF/F), CPT and CP violation (Rez, Imz), 
and CP and T violation (ImAcp, I'z/pI) arc extracted 
from a single fit of both fully reconstructed CP and flavor 
events, tagged and untagged. This provides the sensitiv- 
ity required to separate all effects we seek from asymme- 
tries in detector response and from potentially obscuring 
correlations in the decays of the two B mesons. The 
results are 

sgn(ReAcp)AF/r = 

-0.008 ± 0.037(stat.) ± 0.018(syst.)[-0.084, 0.068] , 

Wp\ = 

1.029 ± 0.013(stat.) ± 0.011(syst.)[ 1.001, 1.057] , 
(ReAcp/|Aa.|)Rcz = 

0.014 ± 0.035(stat.) ±0.034(syst.)[-0.072, 0.101] , 
Imz = 

0.038 ± 0.029(stat.) ± 0.025(syst.)[-0.028, 0.104] . 

The values in square brackets indicate the 90% 
confidence-level intervals. When estimating the lim- 
its we also evaluate multiplicative contributions to the 
systematic error, adding them in quadrature with the 
additive systematic uncertainties. Fig. 4 shows the 
results in the (|g/p|-l, jzj) plane, compared to the 
BABAR measurement of \q/p\ made with dilepton events, 
\q/p\ = 0.998 ±0.006 ±0.007 [11], and to the Standard 
Model expectations. The region shown for this analysis 
is obtained by simulating a large number of experiments 
using the measured covariance matrix for the parameters 
Rez, Imz, and \q/p\, and is constrained to lie within the 
physical region |z| > 0. The three-dimensional distribu- 
tion in Rez, Imz, and \q/p\ is projected onto the two di- 
mensions |z|^ and \q/p\. The boundary is then chosen to 
exclude the maximal region. For simplicity in the figure, 
we display |z| rather than |z|^. The dilepton measurement 
constrains \q/p\ without assumptions on the value of |z|. 
The region in this case is obtained from the Ax^ = 1 



limits for this single variable. 
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FIG. 4: Favored regions at 68% confidence level in the 
(jg/pj — 1, |z|) plane determined by this analysis and by the 
BABAR measurement of the dilepton asymmetry [11]. The 
axis labels reflect the requirements that both CP and T be 
violated if 7^ 1 and that both CP and CPT be violated if 
|z| 7^ 0. The Standard Model expectation for \q/p\ is obtained 
from Refs. [8, 12, 13]. 

Assuming CPT invariance the results are 

sgn(ReAcp)AF/F = 

-0.009 ± 0.037(stat.) ± 0.018(syst.)[-0.085, 0.067] , 

Wp\ = 

1.029 ± 0.013(stat.) ± 0.011(syst.)[ 1.001, 1.057] . 

These results can be used to set constraints on the com- 
plex ratio F12/M12 when CPT invariance is assumed, as 
shown in Fig. 5. Ellipses in the upper figure enclose the 
favored regions determined from the sgn(Re Acp)AF/r 
and \q/p\ measurements of this analysis with z fixed to 
zero. Solid contours show the results assuming Re Xcp > 
(as expected in the Standard Model based on other 
experimental constraints), while dashed contours are for 
Re Acp < 0. Inner (outer) contours represent 68% (90%) 
confidence-level regions for two degrees of freedom. The 
lower figure is an enlargement of the region around the 
origin of the complex ri2/Mi2 plane. The black region 
close to the origin of the complex plane in the upper and 
lower figures shows the predictions of Standard Model 
calculations when all available experimental inputs are 
used to constrain the ratio of CKM matrix elements 
i^cb^cd) / i^tb^td) ■ The bands in the lower figure are cal- 
culated using only the constraint sin 2/3 = 0.741 ± 0.075 
obtained from the BABAR measurement with CP eigen- 
states like J/ipK'^ [5]. 

The decay-rate difference results can alternatively be 
expressed normalized to the mass difference Am. Using 
the world-average value of Am [31], the result allowing 
for CPT violation (z free) is 

sgn(ReAcp)AF/Am = 

-0.011 ± 0.049(stat.) ± 0.024(syst.)[-0.112, 0.091] , 
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FIG. 5: Constraints at 68% and 90% confidence level on 
the complex ratio T12/M12 of the effective Hamiltonian off- 
diagonal matrix elements governing neutral- B-meson oscil- 
lations as determined from the sgn(Re Acp) AF/F and \q/p\ 
measurements of this analysis with z fixed to zero, compared 
to predictions of Standard Model calculations when other ex- 
perimental inputs are used. The lower figure is an enlarge- 
ment of the region around the origin. The bands in the 
lower figure are calculated using only the constraint obtained 
from the BABAR sin 2/3 measurement with CP eigenstates like 
Jf-tpK*^ [5]. The fading out of the bands away from the ori- 
gin indicates that these predictions are only valid for small 
IFia/Mial- 



and with CPT invariance (z = 0) 

sgn(ReAcp)Ar/Am = 

-0.012 ±0.049(stat.)± 0.024(syst.)[-0.113, 0.090] 



The parameters Am and ImAcp/|Acp| are free in the 
fit, so that recent i?- factory Am results [1-4] and our 
sin 2/3 analysis based on the same data sample [5] pro- 
vide a cross-check. The value of the CP- and T-violating 
parameter lmXcp/\\cp\ increases by -1-0.011 when CPT 
violation is allowed in the fit. This change is equal to 
15% of the statistical uncertainty on Im Acp/|Acp| and is 
consistent with the correlations observed in the fit with 
CPT violation. 

The results are consistent with Standard Model ex- 
pectations and with CPT invariance. To date, these are 
the lowest limits on the difference of decay widths of 
mesons and the strongest test of CPT invariance outside 
the neutral-kaon system [19]. If we express the CPT lim- 
its as ratios of the CPT- violating to the CPT-conserving 
terms we have 

M< 1.0x10-", 

m 

ST 

-0.156 < — < 0.042 

at the 90% confidence-level. The limit on CP and T vio- 
lation in mixing is independent of and consistent with our 
previous measurement based on the analysis of inclusive 
dilepton events [11]. All the other results are also con- 
sistent with previous analyses [4, 9, 10, 20, 31-33]. All 
these measurements were obtained with more restrictive 
assumptions than those used here. While the Standard 
Model predictions for AT and \q/p\ are still well below 
our current limits and no CPT violation is anticipated, 
higher precision measurements may still bring surprises. 
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APPENDIX A: EFFICIENCY ASYMMETRIES 

The use of untagged data is essential for determin- 
ing the asymmetries in the tagging and reconstruction 
efficiencies. To indicate how the various samples enter 
we provide a simple example using only time-integrated 
quantities. In practice we use a time-dependent analysis, 
which gives better precision because it uses more infor- 
mation. 

Suppressing the indices for the tag category index a 
and the signal or background component j, and writing 

the reconstruction efficiencies as p = p^^, p = p^- and 



the tagging efficiencies as r = 
reads 



II = 



P + P 

T — T 

T + T 



r%\ Eq. (41) 



(Al) 



Using the numbers of signal events that are tagged and 
have a reconstructed (X), those tagged and having 

a (F), those untagged with a reconstructed (Z), 
and finally those untagged with a reconstructed B^ {W), 



we can determine the required asymmetries [29]. To see 
this, note that if the total number of B'-'B^ pairs is N, 
and neglecting AF, \q/p\ — 1, and z corrections, there are 



Nu = iV(l + [l/(l + x2)])/2 
unmixed events (i.e., B^B'^) and 

= 7V(l-[l/(l + x2)])/2 



(A2) 



(A3) 



mixed events (i.e., or B^B°), where x = Am/F. 

Then we have 



X = pTNm/2 + pfNu/2 , 

Y = pTNra/2 + pTN^/2 , 

Z = p(l - t)7V„/2 + p(l - T)iV„/2 

W = p(l-T)7V^/2 + p(l-r)iV„/2 



Setting U ^ X + Z &ndV = Y + W,wc&nd 



u + v 



[Y/V) + {X/U) 



(A4) 



(A5) 



Corrections to these equations have to be applied due to 
non-zero values of AF, \q/p\ — 1 and z. The use of un- 
tagged events is therefore essential to the determination 
of V and p. 
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